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Abstract. We use the BH masses deduced from the empirical relation of Kaspi et al. (2000) between the size 
of the Broad Line Region (BLR) of Active Galactic Nuclei (AGN) and the optical luminosity, to compute their 
accretion rate in four samples of AGN, assuming that the optical luminosity is provided by the accretion disc. We 
QQ ' show that Narrow Line Seyfert Galaxies 1 (NLSls) accrete at super-Eddington rates, while their luminosity stays 

of the order of the Eddington limit. We take into account the possibility of a non-viscous energy release inversely 
proportional to the square of the distance in the gravitationally unstable region of the disc emitting a fraction of 
^ I the optical luminosity. It leads to a smaller accretion rate and to a redder continuum than a standard disc, which 

agrees better with the observations. The observed bolometric luminosities appear to saturate at a few times the 
Eddington luminosity for super-Eddington accretion rates, as predicted by slim disc models. They favor a Kerr 
BIf rather than a Schwarzshild one. Even when the accretion rate is super-Eddington, it stays always of the order 
of a few M0 /yr, irrespective of the BH mass, indicating that the growing of the BH is mass supply limited and 
' therefore regulated by an exterior mechanism, and not Eddington limited. The mass of the BH increases by one 

, order of magnitude in a few 10^ years, a time smaller than that necessary for changing the bulge mass. This is in 

agreement with recent claims that the BHs of NLSls do not follow the same black hole - bulge relation as other 
galaxies. Since they represent about 10% of AGN up to a redshift of 0.5, these "super-active" phases should play 
an important role in shaping the mass function of local BHs. We finally discuss the possibility that the masses 
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Q ' could be systematically underestimated due to an inclination effect, and we conclude that it could indeed be the 

^ , case, and that the accretion rates could thus be strongly overestimated in a small proportion of objects, possibly 
C/3 ■ explaining the existence of apparently extremely high accretors. 

a 
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1. Introduction and rationale related, so it is generally assumed that their co-evolution is 

mainly the result of merger events within the hierarchical 

The evolution of massive black holes (BHs) in relation • n j. j. r j.- /tt 1 at i 

^ ' scenario of large structure formation (Haehnelt, Natarajan 

with their host galaxy is presently intensively debated. ^ ^^^^ ^ggg^ Kaufman & Heahnelt 2000, Menou, Haiman 

Massive black holes seem present in all galactic nuclei, in- ^ Narayaman 2001, Hatziminaoglou et al. 2003). 

dependently of their level of activity. In about 40 inactive , . ... , ... 

1 1 • ii • r J 1 J. However this scenario begins to be questioned sen- 
nearby galaxies, their mass was found proportional to the , r ■ ■ 1 T Tcr ^ i ■ i n titt 
1 • r it- u 1 -c V, J. 1 r-Mt ■ j. ously. It IS indeed difhcult to explain how smaller BHs 
luminosity of the bulge of the host galaxy (Magorrian et , , , , ^ . _ , 

1 1noo^ p i^/r -j-i /r,ooo\ J /-i ui, ji 4- i grow at lower redshifts and more massive ones at higher 

al. 1988). Ferrarese & Merritt (2000) and Gebhardt et al. ,,.„ „ . , , , , 

/onnn ^ i. J j.u x uj. w • j. x rcdshift. So Marcom et al. (2004) propose that local BHs 

(2000a) showed that a tighter relation exists between the . , , . , • ^ , i • / a ^^tn i 

c 4.1. otT ji^ J 4.1, J- • 1-4. f 4-v, grow mamly during Active Galactic Nuclei (ACjJN) phases, 

mass of the BH, M , and the dispersion velocity ctb of the . , , . , , t^tt • i i 

, , rni, 1 r 1 4-- • 4.-11 J 1, 4. J J 4.U 1 his raises immediately the question whether BHs m local 

bulge. Ihe slope of the relation is still debated, and the , , . „ „ , t^tt ;i , i • 

, r m • X 1 rr^nnn\ ■ 1 1 AGN aud lu quasars follow the same BH/bulge relation- 
recent work of Iremame et al. (2003) gives a value close , . i , • / o 
, . „ 1 1 • i- r • 1 4-- 1 ship as other galaxies, 
to 4. beveral mechanisms accounting for this relation have 

been proposed (Silk & Rees 1998, Umemura 2001, King The BH masses in AGN are not determined like in in- 

2003). When ctb is expressed in terms of the bulge mass, active galaxies by the study of the stellar rotation curve 

it leads to M ~ 0.002M(Bulge). It is thus clear that the close to the center. In about 40 AGN, they are deter- 

growth of the BH and the evolution of the host galaxy are i^ii^^d directly through reverberation mapping (Wandel 

et al. 1999, Kaspi et al. 2000), which yields an empiri- 
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Broad Emission Line Region (BLR), and then to the BH 
mass, using the Full Width at Half Maximum (FWHM) 
of the broad lines as a surrogate of their dispersion veloc- 
ity and assuming that the BLR is gravitationally bound 
to the BH, an assumption confirmed by detailed studies 
(Peterson & Wandel 1999 and 2000). In the other AGN the 
BH masses are determined indirectly assuming that the 
same relations hold. Wandel (1999) showed that Seyfert 
galaxies have lower BH to bulge mass ratios than inac- 
tive galaxies, but the revision of the Magorrian relation 
leads to conclude finally that it is not the case (Laor 2001, 
Wandel 2002, Gcbhardt et al. 2000b). 

However the status of Narrow Line Seyfert 1 galax- 
ies (NLSls) is not well established in this context. NLSls 
constitute about 10% of Seyfert nuclei and quasars up 
to a redshift of 0.5 (Wilhams, Pogge, & Mathur 2002). 
Though they are known since a long time (Osterbrock 
& Pogge 1985), their nature is still not well understood. 
Besides the "narrowness" of their broad lines, these galax- 
ies share common properties, such as strong FeH permit- 
ted lines and weak forbidden [OIII] lines, a strong X-ray 
variability and a big soft X-ray hump (see several reviews 
in Boiler et al. 2000). Mathur, Kuraszkiewicz & Czerny 
(2001) suggested that the BH/bulge mass ratio is smaller 
in NLSls, and Wandel (2002) found that M ~ lO'^ to 
10~^M(Bulge), a smaller value than for broad line AGN 
(BLSls). Both papers are based on a very limited sam- 
ple, and are prone to statistical uncertainties. Moreover, 
in NLSls the bulge mass is generally not deduced from 
the stellar dispersion velocity but from the width of the 
[OIII] 5009 line assumed to be proportional to it, follow- 
ing a suggestion of Nelson and Whittle (1996) for Seyfert 
1 galaxies (actually Wandel (2002) used direct measure- 
ments of the bulge luminosity). Wang and Lu (2001) ar- 
gued that the [OIII] width is not accurately determined in 
NLSls, owing to the weakness of the line and to the pres- 
ence of a blue wing, both effects leading to overestimate 
(t([OIII]) and therefore the bulge mass. However Grupc 
& Mathur (2003) confirmed the previous result of Mathur 
et al. (2001) with a complete X-ray selected sample of 
NLSls, even when taking into account the presence of the 
blue wing of the [OIII] line, and she claims that NLSls 
occupy distinct regions in the BH/bulge mass relation. 
Botte et al. (2004) do not confirm this result, and from a 
study of the photometric properties of the host galaxies 
they find that the NLSl galaxies seem to share the same 
BH/bulge mass relation as ordinary Seyfert, and simply 
occupy the lower ranges of the M — M (Bulge) plane. Bian 
and Zhao (2003) came to an opposite conclusion, based 
also on the bulge luminosity (we recall that the relation 
deduced from the bulge luminosity and the host properties 
is more dispersed than that deduced from the dispersion 
velocity), but found that NLSls do not follow the ordinary 
relation when using the [OIII] line as an indicator of the 
dispersion velocity (Bian fc Zhao 2004). Finally Botte et 
al. (2004) show that there is a smooth relation between 
the BH mass vs. the bulge luminosity for different classes 
of AGN, while there is a jump between the BH mass v.s. 



the [O HI] width. The latter finding is consistent with 
what was claimed by Grupe & Mathur (03) and by Bian 

& Zhao (04). 

One sees that the problem of the BH/bulge mass rela- 
tion in NLSls is presently highly controversial. It has im- 
portant cosmological consequences. If BHs in NLSls arc 
undermassive with respect to their host bulge, it would 
imply that these galaxies are "young", in the sense that 
they are still in the process of building their BH. It would 
mean that BHs and galaxies do not evolve concomitantly 
(Mathur 2000, and Grupe & Mathur 2003). We will show 
here that there is a strong reason to believe this is true, 
because NLSls seem to be accreting at super-Eddington 
rates and therefore the time scale for the growing of their 
central black holes could be extremely short. 

It is now widely admitted that NLSls are radiating 

close to the Eddington luminosity LEdd- This result is 
simply obtained from the mass-luminosity-FWHM rela- 
tions mentioned above. A few objects might even have 
super-Eddington bolometric luminosity, depending on the 
conversion factor used to transform the optical-UV lumi- 
nosity into a bolometric one, and on the adopted Hubble 
constant, but it never exceeds a few ^Edd- From this re- 
sult many people assuming that the efficiency factor for 
conversion of mass into energy is constant and of the order 
of 0.1 deduce that these objects are also accreting close to 
their Eddington limit. 

But why would it have to be so? Super-Eddington ac- 
cretion is indeed theoretically allowed. Near the BH, the 
gas forms an accretion disc, which is supposed to emit the 
"Big Blue Bump" (BBB). The accretion rate and the BH 
mass determine the spectral distribution and the fiux of 
the BBB. It is thus possible to determine the accretion 
rate when the mass is known. It was performed by Collin 
et al. (2002, hereafter referred as C02), using the sam- 
ple of Kaspi et al. (2000) for which the BH masses are 
deduced from reverberation mapping, and assuming that 
the optical luminosity is provided by a standard accretion 
disc (once the luminosity of the underlying galaxy has 
been subtracted). They found that a fraction of objects is 
accreting at siiper-Eddington rates, while their optical lu- 
minosity stays lower than or of the order of the Eddington 
luminosity. Actually, when the accretion rate is close to, 
or larger than the Eddington limit, accretion close to the 
BH does not proceed through a "thin" , but a "slim" disc 
whose cooling time is larger than the viscous time, so en- 
ergy is advected towards the BH before being radiated. 
The mass-energy conversion eflflciency 77 thus decreases as 
the accretion rate increases, and the luminosity increases 
only logarithmically with the accretion rate (Abramowicz 
et al. 1988, Wang et al. 1999, Fukue 2000, Mineshige et 
al. 2000, Wang & Netzer 2003, Kawaguchi 2003). The 
emission of such a disc is characterized by a soft X-ray 
bump as those observed in NLSls. Kawaguchi (2003), 
and Kawaguchi, Pierens & Hure (2004, hereafter called 
KPH) have confirmed that the overall Spectral Energy 
Distribution (SED) of the two most super-Eddington ac- 
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cretors are well fitted by the emission of a slim disc. 
Finally, Wang (2003) noted that super-Eddington accre- 
tion should lead to a limit relation between the BH mass 
and the FWHM of the lines, and he found several ob- 
jects satisfying this relation, indicating that they radiate 
close to their Eddington luminosity, but accrete above the 
Eddington limit. 

There were only a few NLSls in the Kaspi et al. sam- 
ple studied in C02. Moreover the sample is not statistically 
complete since half of the objects are nearby Seyfert nuclei 
chosen mainly for their high degree of variability. The re- 
cent release of several complete samples including a large 
number of NLSls, and the renewed interest for these ob- 
jects since a few years, motivated us to conduct the same 
study on these new samples. While only standard discs 
were assumed in C02, here we take into account the devi- 
ation from the standard disc due to the disc self-gravity, 
which is particularly important in super-Eddington ob- 
jects (cf. KPH). We use also the slim disc model to com- 
pute the bolometric luminosity as a function of the ac- 
cretion rate. We finally discuss some observational conse- 
quences not envisioned in C02. The model can account for 
the fact that the optical- UV continuum of NLSls is red- 
der than that of ordinary Seyferts (Constantin & Shields 
2003). The variation of the bolometric luminosity with 
the accretion rate agrees with the slim disc model. It ex- 
plains why the FWHMs of the broad lines are larger than 
700km/s. 

In this paper, we only want to show some general 
trends and draw qualitative conclusions concerning the 
accretion rates of NLSls, using rough theoretical models 
of accretion discs and applying them to entire samples. 

Finally we insist on the fact that all along this paper 
we accept the commonly admitted statement that the nar- 
rowness of the lines of NLSls is not due to an inclination 
effect, i.e. that NLSls do not constitute a sample of normal 
Seyfert 1 nuclei whose broad line region is a rotating disc 
seen almost face-on. In this case, it is clear that the masses 
derived from the reverberation mapping formulae would 
be strongly underestimated, and consequently their lumi- 
nosity (in terms of Eddington luminosity) overestimated. 

In the following section, we recall first how BH masses 
are determined and we present the samples. We discuss 
the explanation of the empirical relation between the lu- 
minosity and the size of the BLR. In Section 3, we sum- 
marize the theoretical model. Section 4 is devoted to a 
discussion of the results, and in the last section we discuss 
the alternate possibility that the masses of NLSls could 
be underestimated and the accretion rates overestimated. 



2. Determination of the BH masses 

2.1. The empirics! mass-luminosity relation 

Reverberation mapping studies allowed to determine the 
size of the BLR in about 40 objects. It lead to the dis- 
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Fig. 1. Respectively L(5100) (top) and L(5100)/LEdd 
(bottom) versus the FWHM for all samples. The black 
(respt. red) triangles indicate the objects with L(5100) 
< 0.5 lO*"* ergs/sec for all samples (respt. the Kaspi et al. 
sample). The object lying much below the others is NGC 
4051. 

covery of a correlation between the radius of the re- 
gion emitting the H/3 line, which we will call i?(BLR), 
and the monochromatic luminosity at 5100A, L(5100) = 
i^L,,(5100) (Kaspi et al. 2000): 



i?(BLR) = 32.9 X L(5100)°4'^ It days 



(1) 



where L(5 1 00)44 is expressed in lO"'^ erg/s- Though there 
is some uncertainty in the functional form of the rela- 
tion (cf. Laor 2003, Netzer 2003), all recent papers adopt 
this relation to compute i?(BLR) in quasars and Seyfert 
galaxies, when it has not been determined by reverbera- 
tion mapping. 

It is now well demonstrated that the broad H/3 emit- 
ting region is gravitationally bound to the BH (Peterson 
& Wandel 2000). This gives another relation, Mbh = 
R(BLK)V^/G, where G is the gravitational constant. V 
is generally taken equal to VS/2 FWHM, corresponding 
to BLR clouds in random orbital motion. The relation 
becomes, using Eq. ^ 

Mbh = 5.8 10^ x i?(BLR/lt days) x (FWHM) 2000 Mq, (2) 

where (FWHM) 2000 is the FWHM of the H/3 line expressed 
in 2000 km/s (we choose this value since NLSls are de- 



4 



Suzy Collin and Toshihiro Kawaguchi: The accretion rate in NLSls 



fined by FWHM< 2000 km/s). Using Eqs. [T] and [3 one 
gets a relation between AIbh and L(5100) which allows 
to determine M^h as a function of the optical luminosity 
and the FWHM, without the need to know the size of the 
BLR. We stress however that the use of the FWHM as a 
surrogate of the dispersion velocity can lead to a system- 
atic underestimation of the mass, if the BLR is a relatively 
flattened structure dominated by rotation, in which case 
the inclination of the system would play an important role 
(see Section 5). 

These relations have important consequences. If one 
assumes that Lboi ^ 10 L(5100), a canonical value for the 
quasar continuum (cf. Elvis 1994, Laor et al. 1997), one 
gets from Eqs.nand[21 

i?Edd- 0.35L(5100)°4'(FWHM)-'oo (3) 
= 0.28M7"-43(p^jjM)^(,of It days, 

where we call i?Edd the Eddington ratio, i.e. the ratio of 
the bolometric upon the Eddington luminosity LEdd = 
1.5 10**^M7, and Mj the BH mass expressed in lO'^ M©. 
It is obvious from this relation that NLSls have larger 
Eddington ratios than BLSls for a given BH mass. 

2.2. Comments on the luminosity-size relation 

There are several possible explanations for this relation. 
Line emission can be suppressed by dust beyond the radius 
of sublimation, which corresponds to a given heating flux 
oc Lboi/-R^ (Netzer & Laor 1993). But this constraint pro- 
vides only an outer boundary of the BLR. Nicastro (2000) 
proposed that clouds are formed in a wind above the disc, 
close to the transition region between the gas and the ra- 
diation pressure dominated zones of the disc. However the 
size of the BLR depends both on the BH mass and on the 
luminosity, while the observations give only a luminos- 
ity dependence. The striking similarity of AGN spectra 
led also to the idea that the "ionization parameter" (i.e. 
the radiation pressure to gas pressure ratio or the photon 
density to gas density ratio, cx L\,o\/{nE?), n being the 
electron number density) is constant among all objects. 
Actually the size-luminosity relation rather implies that 
the product of the density with the ionization parameter 
is constant. This is consistent with the so-called "LOG" 
model. 

In 1995, Baldwin et al. proposed that the observed 
spectrum of AGN is simply a consequence of the ability of 
a photoionized medium to reprocess the underlying con- 
tinuum "as long as there are enough clouds at the correct 
radius and with the correct gas density to efficiently form 
a given line" . In this "Locally Optimally Emitting Glouds" 
(or LOG) model, each line is emitted preferentially at an 
appropriate ionizing flux L/4ttR^ corresponding to a given 
distance from the source ^. According to the grid of pho- 

^ This is actually closely related with the old idea of line 
saturation due to thermal quenching (Ferland & Rees 1988, 
CoUin-Souffrin & Dumont 1989) 



toionized models pubhshed by Korista et al. (1997) the 
"optimal" ionizing flux Foptimai for the the H/3 line does 
almost not depend on the density and on the spectral dis- 
tribution of the ionizing continuum. It is of the order of 10* 
erg s~^ cm^^. This means that as long as there are clouds 
in a large range of radius with the appropriate density (i.e. 
between 10^ and 10^'* cm~'^) the ionizing continuum will 
be reprocessed in the H/3 line with a maximum efficiency 
at an optimal distance -Roptimai ~ 2 lO^^LP^^ 44 cm, where 
Lion is the ionizing luminosity. From the Gru03 sample, 
one gets L(5100)~ 0.1 Lgjj (precisely L(5100)= 0.21 L^-^^ 
and Lboi = 17 L(5100)^'^^, with a correlation factor of 
0.9). Thus the observed relation transforms into i?BLR ~ 
a few lO^^L?;^ 44 cm, which is similar to the relation ex- 
pected for the LOG model (the ionizing luminosity being 
slightly smaller than the bolometric luminosity). 

So the only necessary condition for the observed rela- 
tionship is the existence of clouds within a broad range 
of density at a radius smaller than the typical distance of 
the BLR, say lO-^i^c- Collin & Hure (2001) suggested that 
such clouds form above the gravitationally unstable region 
of the disc. Since the disc becomes gravitationally unsta- 
ble at small radii compared with the size of the BLR (cf. 
later), this condition is satisfied. The BLR clouds would 
thus constitute simply the outer part of the region emit- 
ting the optical continuum. Laor (2003) objected to this 
idea that "since all accretion discs must become gravita- 
tionally unstable far enough from the center, this mecha- 
nism does not provide a natural explanation for the appar- 
ent absence of a BLR in some Agn" . But there are actu- 
ally several possible explanations for the absence of BLR. 
For instance, in low luminosity objects, it can be due to 
the suppression of the ionizing radiation in an Advection 
Dominated Accretion Flow (ADAF). It can also simply 
be caused by the absence of adequate physical conditions 
in the gravitationally unstable disc, like a too high or to 
small density. 

2.3. The samples 

We use two complete samples including both NLSls and 
BLSls. 

The recent data release of the Sloan Digital Sky Survey 
(SDSS) allowed Boroson to build an homogeneous sam- 
ple of 107 low-redshift radio-quiet QSOs and Seyfert 1 
galaxies (Boroson 2003). It is aimed at comparing the 
BH masses determined from the empirical relations with 
those deduced from the dispersion velocity of the [OIII] 
line, used as a surrogate of the stellar velocity dispersion. 
About one third of objects are NLSls in this sample. It 
allows to study a large range of masses and luminosities. 
For each object the redshift, the FWHM(H/3) and the BH 
mass are given, and we deduce the optical luminosity at 
emission from Eqs.QandlSl We call this sample Bor03. 

The second one is a complete sample of X-ray selected 
AGN (Grupe et al. 2003). According to the selection pro- 
cedure, about half of the objects are NLSls. L(5100) is 
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Fig. 2. These figures display m as a function of the BH masses for the four samples. The squares give m computed 
according both to the standard disc model (open squares), and to the self-gravitating disc model with a viscosity 
parameter a = 0.01 (filled squares). The crosses give L(5100)/Lodd, and the crosses with open circles mark the 
NLSls. The two thick solid lines delineate the position of m for the NLSls. The two horizontal lines correspond to 
M/MEdd = 1, where Msdd = -^Edd/ (^c^); in the case of a Schwarzschild BH and of a Kerr BH. The black (respt. red) 
triangles indicate the objects with L(5100) < 0.5 lO*"* ergs/sec (respt. for the Kaspi et al. sample). Note that the VVG 
sample consists only of NLSls, thus circle symbols are not shown. 



given, but for an empty universe, so we made the con- 
version to go = 0.5. We call this sample Gru03. It is 
particularly interesting for us as it gives an estimate of 
the bolometric luminosity of the objects based on the ob- 
served spectral energy distributions, which we will be able 
to compare with our models. 

We use also two other heterogeneous samples. Wang 
& Lu (2001) deduced L(5100) from the B-magnitude us- 
ing the Veron-Cetty et al. (2001) sample, which contains 
59 NLSls, and they estimated the BH masses using the 
previous empirical relations. After rejection of a few ob- 
jects for which the FWHM are controversial, the sample 
was reduced to 54 NLSls. We call it the VVG sample. We 
also used an heterogeneous sample of soft X-ray selected 
AGN (Grupe et al. 1998, 1999), which has the advantage 



of giving optical indices useful to check our models. We 
also made the conversion from go = to go = 0.5. We call 
it Gru99. Note that a few objects are also in Gru03. 

The samples have not been corrected for the stellar 
contribution of the host galaxy to the optical luminos- 
ity. It is certainly important for low luminosity AGN, but 
not when the optical luminosity is larger than a few 10^^ 
ergs/sec. In the following we will distinguish or suppress 
all these weak objects from the samples, so we can be fairly 
confident that the results will not be contaminated by the 
host galaxy. 

Fig. n] displays respectively L(5100) (top) and 
L(5100)/LEdd (bottom) versus the FWHM for all samples. 
We note immediately the strong difference between these 
two graphs. While the first one shows a very loose cor- 
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Fig. 3. Rsg/ Rg versus M for all samples, for the self- 
gravitating disc with a = 0.01. The black triangles in- 
dicate the objects with L(5100) < 0.5 10'*'' ergs/sec for all 
samples. 

relation, corresponding to the absence of low luminosity 
objects with large FWHMs and of high luminosity objects 
with small FWHMs, the second one shows a tight correla- 
tion with a slope equal to -2, which is expected according 
to the first line of Eq. 3. The black triangles indicate the 
objects with L(5100) < 0.5 lO''^ ergs/sec: note that these 
low luminosity objects share the same relation as the oth- 
ers. 

In an aim of comparison, we have added on these fig- 
ures the objects where the BH masses have been deter- 
mined directly by reverberation mapping (we call these 
objects the "Kaspi et al. sample", though half of them 
were not observed by Kaspi et al. 2000). They span the 
same range of luminosities as the other samples. But first, 
they show a looser correlation between L(5100)/LEdd and 
the FWHMs; it is expected as the determination of the 
mass in the other objects makes use of an exact relation 
L — i?(BLR), not taking into account its error bars. And 
second, the relation should be extrapolated to values of 
the mass and of the Eddington ratio smaller by a factor 
of 5. It should be kept in mind in the following analysis. 
Note that the values of the luminosities used in this figure 
correspond to Hq = 75 km/sec/Mpc, while C02 assumed 
Ho = 50 km/sec/Mpc. 

3. The accretion disc model 

Since more than fifteen years it is widely admitted that 
the "infrared bump" at a few microns and the "Big Blue 
Bump" observed in radio quiet quasars and Seyfert nuclei 
are both due to thermal emission, respectively by hot dust 



heated by the UV-X continuum, and by the accretion disc 
(Sanders et al. 1989). In this picture, the observed "dip" 
at ~5000A in the log(^i^i,) versus logz/ curve corresponds 
to the junction between these two processes, the hot dust 
close to the sublimation temperature (1700K) being un- 
able to radiate appreciably below l/um. In particular the 
idea of an underlying non-thermal power law continuum 
which was invoked in the past and used to model the in- 
frared to UV emission of AGN has been completely left 
over. So the emission at 5000 A should be due entirely to 
the accretion disc, unless another medium can give rise to 
a smooth featureless optical continuum. The problem was 
discussed in C02, and they showed that it would require 
the existence of a very dense, optically thick and relatively 
cold medium. It is difficult to find for such a medium an- 
other location than an optically thick accretion disc. 

For a "standard" thin Keplerian disc where gravita- 
tional energy is released locally through turbulent viscos- 
ity, the effective temperature Toff at a distance R from a 
BH of mass M is: 



(4) 



4 _ 3GMM 

where the non-dimensional factor f{R) takes into account 
the boundary conditions, and is equal to unity at large 
radii (cf. for instance the book of Frank, King & Raine 
2002). 

Each spectral band is emitted around a given radius, 
and the optical band corresponds to a large distance from 
the black hole, typically 10^i?G (-^G being the gravita- 
tional radius GM/c^). At such large radii, the disc is 
dense, relatively cold and optically thick, and its local 
emission spectrum is close to a black body at the tem- 
perature Tcff (cf. Collin 2001; note that it is not the case 
at smaller radii, i.e. in the EUV band). Integrating over 
the disc the Planck law with T cx i?"^/'', one finds for the 
luminosity at a frequency v. 



Rout 



RdR 



.4/3 



Rin 



exp{hv / kT) — 1 



(5) 



where Rin (respt. Rout) is the inner (respt. the outer) 
radius of the accretion disc. 

So it is possible, using Eqs. ^ and to deduce the 
accretion rate when the mass is known. One sees also 
from these equations that L(5100) is approximately pro- 
portional to (MM)^/3 This is not valid for very large 
masses and small accretion rates, where fcr(Rin) is of the 
order of /ifopt, or for truncated discs. 

For super-Eddington accretion rate this picture is 
changed. The radiative efficiency per unit mass accre- 
tion is expected to decrease due to the onset of photon 
trapping (Begelman 1978). As a result, the emergent lu- 
minosity from an accretion flow starts to saturate at a 
few times iEdd (Abramowicz 1988). Self-similar solutions 
with super-Eddington accretion rates (Fukue 2000; Wang 
& Netzer 2003) are only valid inside the photon trapping 
radius, where soft X-ray photons are emitted. However, 
full integration of differential equations from far outside 



Suzy Collin and Toshihiro Kawaguchi: The accretion rate in NLSls 



7 



the photon trapping radius to the vicinity of the central 
BH (Shimura & Manmoto 2003; Kawaguchi 2003) is neces- 
sary in order to discuss the broad-band spectra of NLSls. 
We use the slim disc model for a Schwarzshild BH com- 
puted as in Kawaguchi (2003), which is based on the code 
developed by Matsumoto ct al. (1984). The effect of elec- 
tron scattering (both in opacity and Comptonization) and 
relativistic correction are included. We take the viscosity 
parameter a equal to 0.1. Note that the slim disc is used 
here only to compute the bolometric luminosity. 

Even if the accretion rate is very high (in Eddington 
value) the optical luminosity is still emitted at a large ra- 
dius where the accretion flow is not influenced by advec- 
tion and photon trapping, except in the case of very high 
accretion rates (M > 3 lO^LEdd/c^, cf. KPH), and the 
standard disc model is valid. The only deviation to the lo- 
cal blackbody in the optical region is due to electron scat- 
tering (as modified blackbody, see Czerny & Elvis 1987), 
which distorts the spectrum for super-Eddington accretion 
rates. It is negligible as far as viscosity is small [a < 0.1) 
and the BH mass is small (M < lO'^M©), so the distortion 
is not very important for NLSls (cf. KPH), and we will 
neglect it in this paper. 

However an important fact should not be forgotten, 
which acts also for modest accretion rates but is very im- 
portant for super-Eddington accretion rates. 

At about the distance of the optical emitting region, 
the disc becomes self-gravitating, i.e. the vertical compo- 
nent of the BH gravity becomes smaller than the disc's 
own gravity. This occurs beyond a critical radius Rgg cor- 
responding to a density: 



where Q.k is the keplerian velocity. The disc is then locally 
gravitationally unstable (Goldreich & Lynden-Bell 1965). 
At radii larger than i?sg, the structure of the disc is com- 
pletely unknown. It could break into fragments, which can 
collapse and even form stars, or it can stay at the marginal 
instability limit if it can bo sustained by some extra heat- 
ing mechanism. In all cases the region emitting the 5100 
A flux stays optically thick, and the local black body 
assumption is vaHd (cf. Collin & Hure 1999). 

Koratkar & Blaes (1999) stressed that the standard 
disc model leads to a continuum bluer than the average 
AGN continuum, which has a mean spectral index of 0.3 to 
0.5 (we define the spectral aopt as L^, oc i/~"°p*). Actually 
it is a problem only for small BH masses and large accre- 
tion rates. In the case of large BH masses and small accre- 
tion rates, the optical spcctrum-UV spectrum is emitted 
by the Wien part of the Planck function, and is redder 
than 

Several sources of heating can overcome the gravita- 
tional viscous release beyond the self-gravitational radius. 
The disc can be irradiated by the central source of UV- 
X continuum if it is "flaring" (i.e. if its thickness varies 
more rapidly than the radius). It can be heated by grav- 
itational instabilities (Lodato & Bertin 2003), by the col- 



lisions of clumps (Krolik & Begelman 1988), or by em- 
bedded stars (Collin & Zahn 1999), and cer. In all cases. 
Toff will decrease less rapidly with increasing R than in 
a "standard" disc, and the observed continuum will be 
redder. For instance Soria & Puchnarewicz (2002) fit the 
spectrum of the NLSl 1 RE J1034+396 (this object is 
included in the following computations) by an irradiated 
accretion disc whose scale height to the radius H/ R ratio 
increases rapidly with R, Tea being thus proportional to 
R^^/"^. C02 have shown that in this case, in order to get 
a smooth optical continuum without an intense Balmer 
discontinuity, the density and the optical thickness of the 
irradiated medium should be very large. This is impossi- 
ble with a strongly fiaring disc; a warped thin disc would 
be a more appropriate solution. As we explained previ- 
ously, such a disc would be gravitationally unstable at 
the distance of the region emitting the optical luminosity, 
and most likely very different from a standard one. In the 
case of heating by embedded stars, a very large number of 
massive stars would be necessary to account for the whole 
optical luminosity (Collin & Zahn 1999). 

Since the status of the unstable part of the disc is not 
known, we parametrize these effects by assuming that the 
energy release is proportional to R^^ ^ with /3 smaller than 
3 in the self-gravitating region. In this paper we will as- 
sume the extreme case /? = 2: it corresponds to L^, oc , 
with 7 = 1/2. In the following computations this value is 
used into Eq. 5 instead of Eq. 4 for i? > Rgg, with the 
continuity of the energy release at i?sg- Doing this we ob- 
tain an optical spectral index between -0.3 (corresponding 
to the standard disc) and -1-0.5, depending on the propor- 
tion of the disc which is self-gravitating. It is closer to the 
observed AGN continuum. The effect on the bolometric 
luminosity of this additional energy release is negligible, 
but it increases the computed emission in the optical and 
near- infrared spectral bands, and therefore decreases the 
accretion rate necessary to account for a given optical lu- 
minosity. i?sg is small for small values of a. We have thus 
chosen a relatively small value of the viscosity parame- 
ter (0.01) in order to underestimate Rsg, and therefore 
to underestimate also the accretion rate with respect to a 
standard disc. 

However, we have to take into account the fact that 
the accretion disc cannot extend too much in the self- 
gravitating region, unless a mechanism can act to limit 
the disc density at exactly the marginal instability. Since 
we will see below that the self-gravitation radius is always 
smaller than lO^iJg, we have decided in the following to 
limit the radius of the accretion disc at a value of 10'''i?g. 
It is an arbitrary value, but we have no way to estimate 
the real extension of the accretion disc. Note that the di- 
mension of the BLR is at most of this order in NLSls, and 
it is difficult to accept the idea that the disc extends much 
further out. Note that for such a radius, the gravity of the 
galaxy does not dominate on the BH. 

If the disk is not self-gravitating and extends further 
out, it does not influence the optical emission. Indeed in 
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this case one finds that A(10^i?g) 
which insures that the optical emission is entirely pro- 
duced inside 10^i?g. On the contrary, if the disk extends 
only up to 10^i?g or 10'*i?g, the computed optical emission 
would be smaller than for i?out — W^Rg, and the accretion 
rate would therefore be larger. 

As an accretion disc with a super-Eddington accretion 
rate behaves like a standard disc outside the photon trap- 
ping radius (KPH), we compute Rsg with the same analyt- 
ical approximation as KPH, which gives expressions simi- 
lar to the previous detailed computations of Hure (1998): 



0,5 



R. 



R 



•sg,a 



Rt 



,1/3 



(7) 



'sg,c; 

where i?sg,a, Rsg,h, and i?sg,c are the self-gravitation ra- 
dius in respectively the inner region dominated by radia- 
tion pressure and Thomson opacity, the intermediate re- 
gion, dominated by gas pressure and Thomson opacity, 
the outer region dominated by gas pressure and atomic 
opacity: 
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Fig. 4. The computed optical spectral index aopt mea- 
sured between 4400 and 7000A(rest frame), for the Bor03 
sample, excluding the objects with L(5100)< 0.5 10**^ 
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These expressions depend on the viscosity parameter a. 
We will use a = 0.3, a = 0.1, and a = 0.01. A smaller 
value of a has a more profound influence on the disc struc- 
ture, as it corresponds to a denser standard disc, and 
therefore a smaller value of 

Let us now discuss the consequences of these relations 
in an approximate way. As we shall see later, none of 
the free parameters have a strong influence on the com- 
puted accretion rate, the main quantity that we want 
to determine. We have seen that for a standard disc, 
L(5100) (X (MAf)2/3. Using this relation, and Eqs.Qand 
12 we get: 



puted with the self-gravitating correction, for a equal re- 
spectively to 0.01 (blue circles) and 0.1 (red squares), and 
0.3 (black crosses). We recall that aopt = 0.3 for a stan- 
dard disc. 

cosmological constant, we have made the conversion in 
the aim of uniformity. 

We first draw the attention on a fact which is some- 
times forgotten. Generally it is not the fluxes at Earth but 
the luminosities which are published in the literature, and 
they are computed assuming an isotropic emission. The 
monochromatic luminosity is thus equal to: 



VcL{vc) = AnD'^il + zf X i/c-F(J^c)/Abs(i'o), 



(10) 



rh cx FWHM"''-2^Af 



0,14 



(9) 



where m is the accretion rate expressed in Eddington 
This is actually a very interesting 



M 



units, m = 7 r 

result, which comes from the dependence of the size of 
the BLR on the luminosity and which shows that m de- 
pends almost only on the FWHMs, and very little on the 
BH mass. It is only approximate if the self-gravitating re- 
gion of the disc is large. It means that m can he deduced 
directly from the measurement of the FWHMs alone. 

4. Results and discussion 

We have applied our model to the samples, and we present 
now the results. We use Hq — 75 km/sec/Mpc, and 
go = 0.5. When the luminosities were given for another 



where F is the flux observed at Earth, (respt.i^o) is the 
frequency at emission (respt. at Earth), D is the proper 
distance of the object, z is the redshift, Abs(i'o) the ex- 
ternal (galactic) absorption. But an accretion disc does 
not emit isotropically. The computed monochromatic lu- 
minosity given by Eq. [S] or by its equivalent for the self- 
gravitating region should thus be multiplied by a factor 
2cos(i), where i the inclination of the disc axis on the line 
of sight, to be identified with the published values. 

4.1. Accretion rates 

Figs. 121 display m as a function of the BH mass for the 
four samples, m is computed according both to the stan- 
dard disc model, and to the self-gravitating disc model as 
explained in the previous section. In this latter case, the 
results are shown for a viscosity parameter a = 0.01. In 
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Fig. 5. Comparison between the observed and computed 
optical spectral index ttopt, for the Gru99 sample, exclud- 
ing the objects with L(5100)< 0.5 lO^"* erg/s. The black 
filled circles are the observed values, the open symbols are 
computed with the self-gravitating correction, for a equal 
respectively to 0.01 (blue circles) and 0.1 (red squares). 
We recall that aopt = 0.3 for a standard disc. 

all computations, cos(i) is set equal to 0.75. The objects 
with L(5100) < 0.5 10**'' ergs/sec are indicated on the fig- 
ures. There are only two such objects (actually lying close 
to the limit) in Bor03. The Gru03 sample contains many 
low luminosity objects, but a large number of NLSls are 
above the luminosity limit. 

We see that the self-gravitation correction can decrease 
rh by about a factor three for large values of m, but 
has no influence on small rh. For larger values of a and 
of 7, the difference between the standard and the self- 
gravitating disc would be smaller. So we can consider that 
the two models here correspond to a kind of "error bar" 
on m, for given BH mass and L(5100). Figs. [21shows also 
the "observed" ratio L(5100)/iEdd- We have noted the 
NLSls, and the thick dotted lines delineate the position 
of m for NLSls. NLSls always have BH masses smaller 
than 10* Mq, and they are located in the higher range of 
L(5100)/Ledd and m. It is interesting to note that the four 
samples do not differ except for the range of masses and 
luminosities, though they have been selected quite differ- 
ently. 

Again we added for comparison to these figures the 
results for the Kaspi et al. sample, computed using only 
the standard disc emission (we recall that the results differ 
from C02 because we use here Hq = 75 instead of 50). As 
expected, the extrapolation by a factor 5 in mass range of 
the empirical relationship translates in an extrapolation of 
rh by about a factor 30, as m oc M/Mbh « ^5(00 ^ -^^bh- 



Fig. 6. Accretion rates in Mo/yr for the four samples as 
a function of the BH masses, excluding the objects with 
L(5100)< 0.5 lO**** erg/s, and computed according to the 
self-gravitating disc model with a viscosity parameter a = 
0.01. NLSls are indicated as red dots. 

Several other results appear on these figures. 

First rh increases as the BH mass decreases. On the 
contrary, the ratio L(5100)/LEdd is always smaller than 
0.3, and seems about constant for the NLSls. When ap- 
plying a standard correction Lbol ~ 10 x L(5100), one 
concludes that Lbol saturates at about the Eddington lu- 
minosity, whatever the BH mass. This excludes the ex- 
istence of the large super-Eddington ratios proposed by 
Begelman (2002) due to the photon bubble instability. 
Thus, according to Eq. 3, there should be a lower limit 
to the FWHMs of the order of 1000 M^ '^^ km/s unless 
the empirical relations do not apply to these objects. And 
indeed FWHMs of the order of 100-500 km/s which would 
imply Eddington ratios larger than 10 have never been ob- 
served in Seyfert 1 nuclei. 

Second, the two horizontal lines correspond to 
M/MEdd — 1, where A/Edd = -^Edd/('?c^), in the case 
of a Schwarschild BH (77 ~ 0.057) and of an extremely 
rotating Kerr BH (77 = 0.30). We see that the accretion 
rates of NLSls are always larger than the Eddington rate 
in the case of Kerr BHs, and mostly larger in the case of 
Schwarzschild BHs. 

There are several causes of uncertainties in the results 
(cf. Krolik 2001 and C02), which might introduce errors 
on the BH masses as large as one order of magnitude, 
because one should not forget that even the masses de- 
termined directly with reverberation mapping are known 
with an uncertainty of a factor 3. It seems however im- 
plausible that all the uncertainties would systematically 
act towards an underestimation of the mass and an over- 
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Fig. 7. The Eddington ratio i?edd as a function of m for 
the Gru03 sample, computed with the standard disc (open 
squares) and the self-gravitating disc, a — 0.01 (filled 
squares). The objects with L(5100) < 0.5 10"*^ ergs/sec 
have been suppressed. The two curves correspond to the 
slim disc model, a = 0.1, and respectively a Schwarzshild 
and a kerr BH. 



estimation of the luminosity, avoiding the conclusion of 
super-Eddington accretion rates. Only the uncertainty on 
the correcting factor of the FWHM due to the geome- 
try and kinematics of the BLR could lead to a systematic 
underestimation of the mass, if the BLR is a rotating flat 
structure. It can be large when the objects are seen almost 
face-on. We shall discuss this point in the last section. 

Fig. m displays Rgg/Rc versus M for all samples, for 
the self gravitating disc with a = 0.01. We note that it is 
always quite small (in particular smaller than the BLR, 
which has typical values lO'^ for high BH masses and 10^ 
for NLSls), justifying our previous claim that the BLR is 
always located in, or above, the unstable part of the disc. 
As expected Rsg/Rc decreases with the BH mass, except 
at the high mass limit, and there is a strong correlation 
between the two parameters. 

Although the choice of parameters for the self- 
gravitating disc does not influence strongly m, it has an 
effect on the optical spectral index. As an illustration. Fig. 
0] shows the computed optical spectral index aopt defined 
as Fi, oc zy~"°p' between 4400 and 7000A (rest frame), 
for the Bor03 sample. The computation is performed with 
the self-gravitating correction, for a viscosity parameter 



a equal to 0.01, 0.1, and 0.3. A systematic correction 
E{B — V)—0.05 for the galactic absorption has been ap- 
plied (certainly an underestimation). For a — 0.01, the 
continuum is red except for very broad line objects. The 
trend that broader objects have bluer optical spectra is 
consistent with the observational results of Constantin & 
Shields (2003). The continuum is globally bluer for smaller 
values of a (0.1 and 0.3). We also see that aopt almost 
never reaches the value of the standard disc (-0.33). A de- 
tailed comparison with the observed values is postponed 
to the next paper. 

Fig. |31 shows a comparison between the observed and 
computed spectral indexes for the Gru99 sample, exclud- 
ing the objects with L(5100)< 0.5 10^"* erg/s. According 
to Grupe et al. (1999), the observed values of aopt are 
given with an uncertainty of ±0.4. With a = 0.1, many 
of the computed indices are close to the value of the stan- 
dard disc, while the objects of the samples are particularly 
red, with an average index of 0.8. The agreement is much 
better for the smallest viscosity parameter a = 0.01. The 
very red spectra observed in a fraction of objects might be 
due to intrinsic reddening not taken into account in the 
computed values. If it is the case, it would imply that the 
observed L(5100) is underestimated in these objects, but 
again it is not important for the determination of rh. Note 
that in this sample, NLSls do not seem to have redder 
continua than BLRls. 



It is therefore impossible from this comparison to de- 
cide which are the best values of a and 7 to choose for 
the disc. Our model is clearly oversimplified, and would 
require a more sophisticated parametrization. The only 
conclusion which can be drawn is that a non-standard 
disc with an additional release of energy in its external 
region gives a better fit to the average optical continuum 
of AGN than a standard disc. However, this problem does 
not question the existence of super-Eddington accretion 
rates for NLSls. 

Finally Fig. displays the accretion rates in Mq / yr 
for the four samples, excluding the objects with L(5100)< 
0.5 lO^"' erg/s, and computed according to the self- 
gravitating disc model with a viscosity parameter a — 
0.01. Note that the imposed limit on L(5100) creates the 
sharp limitation on the left side, as AI is proportional to 
7Vf~2/3 ^Jqj. fixed Lopt). The limitation on the right side 
is due to a limitation of m at about 0.03 (perhaps due to 
the fact that the accretion disc changes into an ADAF be- 
low this value). NLSls are indicated as red dots. Despite 
the large values of m of NLSls, we see that the maximum 
accretion rate is of the order of one MQ/yr whatever the 
BH mass. This is a strong indication of an exterior regu- 
lation of the accretion, rather than the self-regulation of 
the disc. Note that it is a modest value when compared 
with the rate of star formation in a starburst nucleus. 
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Fig. 8. Cumulative number of objects (normalized to 
unity) on which an underestimation of the mass by a fac- 
tor smaller than Mreai/Afobs is made, in the conditions 
explained in the text. 

4.2. Comparison with the slim disc model 

It is interesting to compare the observed SED of super- 
Eddington objects with the slim disc model. As we 
mentioned in the introduction, this was done in de- 
tail for the two highest m objects (Kawaguchi 2003; 
KPH; Kawaguchi, Matsumoto, Leighly in preparation; see 
Kawaguchi 2004). and it will be performed for the objects 
of the samples in a future paper. Here we simply com- 
pute the bolometric luminosity, and we compare it with 
the observed values. 

Only the Gru03 sample provides bolometric luminosi- 
ties based on the observed SEDs. Fig. [3 shows the ob- 
served ratio i?Edd versus m for this sample. The low lu- 
minosity nuclei (L(5100)< 0.5 10''^'* ergs/sec) have been 
suppressed. Also shown is the theoretical curves obtained 
for the slim disc model with a Schwarzschild and a Kerr 
BH. These curves depend very little on the BH mass and 
on the viscosity parameter. In spite of the large dispersion 
of the "observations" , it is clear that a majority of points 
lie above the Schwarzschild curve, meaning that the effi- 
ciency of the Schwarzschild BH is insufficient, i.e. a Kerr 
BH with an efficiency of about 0.15 would better fit the ob- 
servations unless there is a systematic underestimation of 
the BH masses. On the other hand, the shape of the curve 
agrees well with the observed points, in particular in the 
"saturation" of i?Edd above m — 10. Three objects reach 
an Eddington ratio of the order of 10, for 50 < m < 1000. 

5. Influence of the inclination on the masses and 
accretion rates 

In all mass determinations, the FWHM is used instead of 
the dispersion velocity. It makes the implicit assumption 
that the velocities are distributed at random in the BLR. 



However, if the BLR is a flat structure dominated by ro- 
tation, the FWHM is proportional to sin(«)VKep, where i 
is the angle between the normale and the line of sight (the 
inclination). It is clear that a small inclination can lead to 
a large underestimation of Vicep and therefore of the mass. 

However the BLR cannot be a geometrically thin disc 
with an exactly Keplerian velocity. Unfortunately its dy- 
namics and its structure are still not well determined 
from detailed reverberation mappings, but we know that 
it should be at least a "thick disc" , with an aspect ratio 
larger than, say, H/R ~ 0.3 {H being the disc thickness 
at the radius i?), since it needs to have a large coverage 
factor of the central source. Such a disc must be sustained 
vertically by a turbulent pressure corresponding to a tur- 
bulent velocity of the order of VKcpH/R. The FWHM is 
then proportional to T^Cep \/ {H/RY + sin(z)2, and the ra- 
tio G between the real mass and the "observed" mass, 
is: 

G = Af,cai/Mobs = ll[{H/Rf + sm{if]. (11) 

We can compute how many objects have a mass un- 
derestimated by a given factor G, assuming that they 
are distributed at random inside an angle iq. We choose 
io = 7r/4 in the following computations, as it is a com- 
monly accepted value for the opening angle of the dusty 
torus in Seyfert 1 (according to the Unified Scheme, cf. 
Antonucci & Miller 1985). The probability of seeing an 
object at an inclination angle i per unit angle interval is 
thus sin(«)/[cos(io) — 1]. The number of objects per unit 
interval of G is: 

m ^ [{H/Rf + s\n{iff 

dG 2[1 - cos(io)]cos(i) ■ ^ ' 

Fig. IHl shows the integral of this expression, i.e. the 
cumulative number of objects (normalized to unity) on 
which an underestimation of the mass by a factor smaller 
than Afroai/A^obs is made, for two values of H/R. Note 
that io plays a non negligible role here, as it contributes 
to increase the proportion of objects having a large G by 
a factor 3 with respect to an isotropic distribution. We 
see that the factor G depends strongly on the aspect ratio 
the BLR. For H/R — 0.3, it can take values as large as 
11, but the number of objects reaching this value is small: 
only about 20% have a G-factor larger than 6, and 1% a 
factor larger than 10. For H/R — 0.5, the maximum value 
of G is only 4, and about 60% have a G-factor of the order 
of or smaller than 2. 

It is important to realize that in this case not only 
nearly face-on objects, but all Seyfert 1 will have their BH 
mass underestimated. This would raise a problem concern- 
ing the BH-bulge mass relationship. 

Would we have thus to modify our conclusions? For 
H/R — 0.3, half of NLSls could have their masses un- 
derestimated by factors of 3 to 10, leading to underesti- 
mations of the accretion rates (in terms of Eddington) by 
factors 10 to 100. It is clearly very important, but still a 
large proportion of objects would be accreting above the 
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Eddington limit, however at a smaller rate. On the other 
hand, it is quite possible that the few extremely high ac- 

crctors arc actually "facc-on" objects, and that their mass 
is indeed underestimated by about one order of magni- 
tude. 

6. Conclusion 

We used the BH masses deduced from the size-luminosity 
relationship to compute their accretion rate in four sam- 
ples of AGN, assuming that the optical luminosity is pro- 
vided by the accretion disc. Thus the empirical relation 
must be extrapolated in a range of masses almost one or- 
der of magnitude smaller than the Kaspi et al. sample. 
We used a simplified disc model, with a paramctrization 
of the energy release in the self-gravitating region to get 
the accretion rate, and the slim disc model in the inner re- 
gions in order to get the bolomctric luminosity. In spite of 
the crudeness of the treatment, this study leads to several 
fairly certain conclusions. 

— NLSls are always accreting at Eddington or super- 
Eddington rates, m can reach 1000, corresponding to 
an accretion rate equal to 60MEdd (for Schwarzschild 
BH) and to 300MEdd (for Kerr BH). 

— Their observed bolometric luminosities "saturate" at 
~ 10 Eddington luminosities, as predicted by slim disc 
models. It explains why there is a lower limit to the 
observed FWHM. 

— The observed value of the bolometric luminosities 
are in better agreement with a Kerr than with a 
Schwarzschild BH. 

— The computed optical spectral indexes agree with the 
observed trend of redder spectra for NLSls than for 
BLSls. 

— And finally the accretion rates have an upper limit of 
about one M0/yr, whatever the BH mass. In partic- 
ular, all NLSls have an accretion rate of this order. 

This is a strong indication for a mass limited supply, 
implying an exterior regulation of the accretion. 

With these results wc are in a position to say now that 
NLSls should have a strong influence on the growth of 
BHs. This is in agreement with the claim by Mathur et al. 
2001, and Grupe & Mathur 2004. Since NLSls constitute 
about 10% of normal Seyfcrt which themselves are about 
2% of inactive galaxies, one deduces that all galaxies spend 
0.2% of their lifetime in the NLSl phase, i.e. 2 10'' years. 
During this time the mass of the BH increases by one order 
of magnitude (Kawaguchi et al. 2004). This could account 
both for the observed large dispersion in the BH/bulge 
mass relation of NLSls, and for the existence of undermas- 
sive BH/bulge ratios during a large fraction of the NLSl 
phase. The increase of the bulge mass could have taken 
place during merger or interaction events. BHs would then 
grow during intense phases of activity after a time delay, 
necessary for accumulating matter in the circumnuclear 
region and for triggering a starburst. In this scenario, the 



overabundance of iron could be easily explained by the 
rapid formation of massive stars and supernovae explo- 
sions in the outer parts of the accretion disc where the 
accretion rate is high (Collin & Zahn 2000, Levin 2003, 
Levin & Belobodorov 2003). The scenario would also ac- 
count naturally for the presence of outflows giving rise to 
the blue wing of the [OIII] line, as super-Eddington accre- 
tion is expected to generate outflows by strong radiation 
fields. 

Though wc have tried to determine a lower limit of the 
accretion rate, two effects can intervene to still reduce it. 
They were both discussed in C02. 

1- the possibility that the accretion rate decreases with the 
radius between the optically emitting region and the BH, 
owing to the creation of a strong outflow due to the radi- 
ation pressure. The accretion rate close to the BH would 
then be just Eddington. In this case, the outflow could 
well be the origin of the [OIII] wing, and could lead to 
the escape of a part of the Narrow Line Region, explain- 
ing the weakness of the [OIII] line. However, one should 
realize that in this case the rate of outflow would have 
to represent 90 or even 99% of the accretion rate, in the 
highest accretors. This seems unrealistic. 

2- The optical luminosity is not provided by the accretion 
disc. Recently King & Pounds (2003) suggested that BHs 
accreting at super Eddington rate produce winds which 
are Thomson thick and can emit a black body spectrum 
providing the Blue Bump of AGN. Pounds et al. (2003) 
indeed report that they have found the signature of such 
an optically thick wind in the X-ray spectrum of the NLSl 
PG1211-I-143. If the existence of such a wind is confirmed 
in other NLSls, then it is clear that the present analysis 
would have to be reconsidered. However let us recall that 
Collin et al. (2002) have shown that very strong conditions 
must be met in such a wind to give rise to the optical-UV 
featureless continuum: it must have both a large density 
(10^^ cm^''), and a Thomson thickness at least of unity. 
Besides, to get the observed luminosity, it should be lo- 
cated far from the center and it should have a large spatial 
extension. It is thus not obvious that the wind observed 
by Pounds et al. (2003) satisfies these requirements. It is 
more likely that its emission is limited only to the EUV 
radiation, and that the optical emission is still due to the 
accretion disc. 

We have assumed all along the paper that the BH 
masses of NLSls are correctly estimated by the empirical 
reverberation relations, even when these relations had to 
be extrapolated by almost one order of magnitude. On the 
other hand, we have accepted the usual assumption that 
the FWHM is a good measure of the velocity in the BLR, 
implicitly assuming that the velocities are distributed at 
random. If on the contrary the BLR is a flat structure 
dominated by rotation, the BH masses of a fraction of ob- 
jects could be underestimated by factors up to one order of 
magnitude and the accretion rates (in terms of Eddington) 
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by two order of magnitudes when they are seen nearly face- 
on. However, since this fraction should be small, we think 
that the scenario described in this paper is qualitatively 
correct. 
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